Preprint typeset using I^-T^X style cmulatcapj v. 5/2/11 



MEASUREMENTS OF STELLAR INCLINATIONS FOR KEPLER PLANET CANDIDATES 

Teruyuki Hirano 1 , Roberto Sanchis-Ojeda 2 , Yoichi Takeda 3 , Norio Narita 3 , Joshua N. Winn 2 , Atsushi 

Taruya 1,4,s , and Yasushi Suto 1 ' 4,6 



ABSTRACT 

We present an investigation of spin-orbit angles for planetary system candidates reported by Kepler. 
By combining the rotational period P s inferred from the flux variation due to starspots and the 
projected rotational velocity Vsvnl s and stellar radius obtained by a high resolution spectroscopy, 
we attempt to estimate the inclination I s of the stellar spin axis with respect to the line-of-sight. For 
transiting planetary systems, in which planetary orbits are edge-on seen from us, the stellar inclination 
/ s can be a useful indicator of a spin-orbit alignment/misalignment. We newly conducted spectroscopic 
observations with Subaru/HDS for 15 KOI systems, whose lightcurves show periodic flux variations. 
After detailed analyses of their lightcurves and spectra, it turned out that some of them are binaries, 
or the flux variations are too coherent to be caused by starspots, and consequently we could constrain 
stellar inclinations I s for eight systems. Among them, KOI-262 and 280 are in good agreement with 
I s = 90° suggesting a spin-orbit alignment, while at least one system, KOI-261, shows a possible 
spin-orbit misalignment. We also obtain a small I s for KOI-1463, but the transiting companion seems 
to be a star rather than a planet. The results for KOI-257, 269, 367, and 974 are ambiguous, and can 
be explained with either misalignments or moderate differential rotation. Since our method can be 
applied to any system having starspots regardless of the planet size, future observations will allow for 
the expansion of the parameter space in which the spin-orbit relations are investigated. 
Subject headings: planets and satellites: general - planets and satellites: formation - stars: rotation 
- techniques: spectroscopic 



1. INTRODUCTION 

The standard formation theory of close-in gas gi- 
ants (hot-Jupiters) suggests that they form outside of 
the so-called "snow-line", located at a few AU away 
from the host star, and subsequently migrate inward 
(e.g iLin et al.l [l99l IChambersI l20ul ILubow fc Idal 
120101) . While migration processes such as disk- 
planet interactions (type I or II migration) predict 
small orbital eccentricities of planets and small stel- 
lar obliquity (i.e. the angle between the stellar 
spin axis and the planetary orbital axis), dynami- 
cal processes including planet-planet scattering and/or 
Kozai cycles might produce large value s for both ec 
centr i city and stellar obliquity (e.g.. IWu fc Murra 1 



20031 iFabrvckv k Tremaindl2007t iNagasawa et al.l l200 
Chatterjee et al.ll2008l iNagasawa fc Idall20liF ~ 



Measurements of the Rossiter-McLaughlin effect (here- 
after, the RM effect), which is an apparent radial ve- 
locity anomaly during a planetary transit, provide in- 
valuable informatio n to better understand the planetary 
migration pr ocess ( Ohta et al.1 120051: IWinn et al.1 120051: 
iNarita et ahl l2007t lAlbrecht et al.l l2007t iHebrard et al.l 
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120081: iTriaud eHII[2010l : IHirano et al.l 1201 lbh . Through 
the RM effect, one can measure the sky-projected angle 
A between the stellar spin axis and planetary orbital axis, 
which is of importance to distinguish among the possible 
migration channels. 

So far, more than 50 transiting systems have been in- 
vestigated to estimate A, and many interesting correla- 
tions among spin-orbit angles and physical properties of 
planets and host stars have been proposed^. I Winn et"al"1 
(|2010al) pointed out that a substantial spin-orbit mis- 
alignment tends to be observed around hot stars (whose 
effective temperatures T c g > 6250 K). This fact might 
be r elated to the tidal evolution of host star's obliqui- 
ties. ITriaudI (|2011| ) also found a correlation between the 
stellar ages and obliquities, claiming that spin-orbit mis- 
alignments are observed around younger systems with 
ages less than 2.5 G yr. This trend is consistent with what 
W inn et aLl (|2010a| ) found, suggesting that tidal interac- 
tions between the host stars and close-in giant planets, 
formed after some dynamical processes such as planet- 
planet scattering, gradually damp the obliquity of host 
stars close to 0°. 

It should be noted that measurements of the RM effect 
are only feasible for rather bright stars (V < 12) with 
giant transiting planets. While detection s of the RM 
effect f or a super-Neptune were reported (jHirano et al.1 
I2011al IWinn et al.l l2010bO . those for smaller planets 
(R p < 0.5Rj) are still challenging. Nevertheless, in or- 
der to discuss planetary formation and migrations, it is 
of great interest to investigate the spin-orbit relations for 
systems with Neptune-sized or even Earth-sized exoplan- 
ets, which are reported to be more abundant than jovian 



7 The list of RM me asurements 
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Fig. 1. — Schematic figure for the spin and orbital axes. In this 
figure, the x — y plane indicates the sky plane and the z axis points 
toward us. The planet's orbital axis is indicated by the blue arrow 
and its projection onto the sky plane is shown by the dashed arrow. 
The red arrow is the stellar spin axis, which is located in the y — z 
plane. 



planets (jBorucki et al.11201 11 [Mayor et al.ll2011l) . 

In order to measure spin-orbit relations for such 
smaller planets, we focus on the Kepler photometry in 
this paper. The Kepler mission is an ambitious and very 
productive project designed to "determine the frequency 
of Earth-size planets in and near the habitable zone of 
solar-type stars". As of March 2012, more than two 
thousand plan etary candidates we r e announced by the 
Kepl er team (jBorucki et al.l 120101 120111: IBatalha et ail 
12012ft . and each of those stars having planetary candi- 
date^) is called Kepler Object of Interest (KOI). Among 
the published Kepler light curves, there are many sys- 
tems that show periodic flux vari ations most likel y due 
to starspots on the stellar surface (jBasri et al.ll201ll ). For 
those systems, a period analysis enables us to infer the ro- 
tational period P s of the star. That rotational rate, along 
with the stellar radius R s , can be directly translated into 
the rotational velocity at the equator of the star, V cq 
in the absence of differential rotation. If we compare 
V cq with the projected rotational velocity Vsinl s esti- 
mated by a spectroscopic observation, we can constrain 
the stellar inclination I s , which is defined as the angle 
between the line-of-sight and the axis of the stellar rota- 
tion (see Figure Q}. Since transiting planetary systems 
nearly have edge-on orbits seen from our location, the 
orbital inclination I Q (the angle between our line-of-sight 
and the planetary orbital axis) should be close to 90°. 
Therefore, a significant deviation of I s from 90° implies 
a spin-orbit misalignment. This method to constrain the 
stellar obliquities will be described in detail in Section [2j 

The statistics of the spin-or bit angles along th e line- 
of-sight has been discussed bv iSchlaufmanl (|2010j) . Us- 
ing the empirical relation among the stellar ages, masses, 
and rotational periods, he statistically discussed the rota- 
tional velocity of the stars hosting transiting planets. He 
found that some of the transiting systems have smaller 
V sin I s than expected for the case of spin-orbit align- 
ment, and pointed out that it is most likely to be evi- 
dence of spin-orbit misalignments along the line-of-sight. 
In contrast to his analysis that relies on an empirical re- 
lation to estimate the rotational period P s of the planet 
hosting stars, we attempt to derive it more directly for 
each of the KOI systems using the precise Kepler pho- 



tometry. In order to constrain I s , we newly conducted 
high resolution spectroscopic observations and obtained 
spectra for 15 KOI systems using the Subaru telescope. 

The remaining sections are organized as follows. In 
Section [5J we describe the basic method to estimate the 
spin-orbit angle along the line-of-sight in detail and dis- 
cuss the pros and cons of the present method in compar- 
ison with the RM effect. We briefly describe the spec- 
troscopic observations with Subaru in Section [3] Section 
2] presents the photometric and spectroscopic analyses 
and their results. The correlations between the stellar 
inclinations and other system parameters are discussed 
in Section [5J Finally, Section [5] is devoted to summary 
and future prospects. 

2. PRINCIPLE 

If there exists a spot on the surface of a star, the 
lightcurve of the star exhibits a periodic variation due to 
the stellar rotation. A period analysis of the lightcurve 
enables us to estimate the rotational period P s of the 
star. Figure [2] shows an example of the Kepler lightcurve 
(KOI-261) and its periodogram. The peak in the pe- 
riodogram most likely reflects the rotational period P s , 
which is estimated as P s = 15.4±0.3 days. Once P s is 
estimated, the stellar inclination I s is estimated by the 
following relation: 



P S (V sin J 5 ) spc 
2ttR, 



(1) 



where (Vsin/ S ) spec and R s are the projected rotational 
velocity and the radius of the star, respectively. Both of 
these quantities arc estimated via spectroscopy under the 
assumption that the star is rigidly rotating. The impact 
of differential rotation will be discussed in Section [5] 

Since the configuration of the transiting system has an 
edge-on orbit seen from us (with the orbital inclination 
I Q > 85°), a small value of I s implies a possible spin- 
orbit misalignment in the system, regardless of the sky- 
projected spin-orbit angle A (see Figure [1]). The 3D angle 
ip between the stellar spin axis and planetary orbital axis 
is associated with J ? , and A by the following equation 
(|Fabrvckv fc Winnl 120091) : 

cos ijj = sin I s cos A sin I a + cos I s cos I Q . (2) 

The present method has several advantages that we 
mention below. Although observations of the RM effect 
have enabled us to discover the spin-orbit misalignment 
for the first time and revealed its possi ble patterns in 
terms of stellar and planetary properties (iHebrard et al. 
20081 INarita et al.l[2009t IWinn et alJl2lMl2010allTr!aud' 



20111 ). the methodology is feasible only for sufficiently 
bright stars (V < 12) with giant transiting planets 
(R p > 0.5Rj). In addition, a measurement of A requires 
spectroscopic observations throughout a complete tran- 
sit. Therefore the observation is time-critical and also 
time-consuming, in particular for those systems in which 
the semi-major axis of the planet is large and the transit 
duration is long. 

In contrast, the current method of measuring spin- 
orbit relations uses the effect of stellar spots induced 
by rotation on the flux of the star, and thus requires 
only one spectroscopic observation independently of the 
size of planets and their semi-major axes. Therefore the 



Measurements of Stellar Inclinations for KOI Systems 



3 



1.005 



ra 1.000 

CD 

s 

£ 0.995 



1.005 



T3 
T3 



1.000 



CD 

a 



0.995 



100 200 300 400 

BJD-2454900 [days] 



500 



4000 F 



3000 



| 2000 

1000 




■4— Jtfc i«fciA 



10 20 30 

Period of rotation [days] 



40 



50 



Fig. 2. — Top. - Quarter by quarter mean normalized raw out- 
of-transit flux of KOI-261 (PDC SAP_FLUX). Includes quarters 
through 6, except quarter 5 where this target was not observed. 
Central.- Detrended flux of the star, obtained with the method 
described in Section [4] Bottom. - Lomb-Scargle periodogram of 
the detrended flux. The solid vertical lines represent the points 
where the power is half of the maximum power. The final value is 
taken to be the mean of these points, and the error to be half of 
their difference. 



present method is more efficient in increasing the num- 
ber of samples, and should expand the parameter space 
in which spin-orbit relations are investigated. 

In systems where spots are present, there is a pos- 
sibility that the planet passes in front of one of them 
during a transit. The recurrence of these spot-crossing 
events at consecutive transi ts proves that the system - 
has a low spin-orbit angle dSanchis-Oieda et all 120111 : 
iNutzman et all 120111 : IDesert et all 120111 ). However, to 
locate the spot anomalies for this spot-crossing method 
requires a high signal-to-noise ratio for the transit 
lightcurve, which is achievable only for large planets or 
very bright stars. Also it would be more suitable for cool 
stars where active regions tend to be larger and produce 
more visible spot anomalies. In addition to this, cool 
stars generally have a lower Vsin/ S , which complicates 
the use of our method. 

On the other hand, measuring the spin-orbit relation 
along the line-of-sight with the present technique has a 
few shortcomings. First, because of the shape of the sine 
function, the uncertainty for I s tends to be larger when 
sin/, is close to unity. Second, it is impossible to dis- 
tinguish the state of I s from that of it — I s , which gives 
exactly the same value of sin/ s . This is a marked con- 
trast with the RM measurements, through which we can 
clearly distinguish between the prograde and retrograde 
orbits at least on the sky. 



Finally, types of stars to which we can apply our 
present te chnique are fairly li mited. According to the 
survey by iRadick et al.l (|1982l) , the flux variations due 
to dark starspots are only seen in Sun-like stars whose 
effective temperatures are up to 6400 K, with a photo- 
metric precision of ~ 3 mmag. Therefore, if we observe 
a periodic flux variation with a large amplitude for an 
early type star hotter than 6500 K, it is more likely to 
be spurious, possibly reflecting a flux contamination by 
a companion star or faint background sources. Therefore 
the present method is complementary to the RM effect, 
in the sense that RM measurements are applicab le to 
even hotter stars (e.g.. iCollier Cameron et al.ll2010D . 

There are some systems in which spots were used to 
estim ate rotational periods of the planet -hosting stars 
(e.g., ILanza et all 120091 IDesert et alJl20Tll for CoRoT-2 
and K epler- 17, respectively). Specifically. IHebrard et all 
(|2011[) applied the present technique to investigate the 
spin-orbit relations to the CoRoT-18 system, and ob- 
tained a weak constraint on the stellar inclination as 
I s = 70° ± 20°. This was in good agreement with 
the small value of the sky-projected spin-orbit angle 
(A = 10° ± 20°), implying a 3D spin-orbit alignment. 
In this paper, we systematically apply this technique to 
some of the KOI systems showing periodic flux varia- 
tions. 

3. OBSERVATION 

We applied the following three criteria to all the KOI 
syste ms in the February 2011 data release ([Borucki et all 
12011 . and selected 15 KOI systems (KOI-4, 42, 100, 
257, 258, 261, 262, 269, 279, 280, 302, 367, 974, 1020, 
1463); 1) the lightcurve exhibits flux variations by a vi- 
sual inspection, 2) the (photometrically estimated) effec- 
tive temperature is higher than ~ 6000 K, and 3) the 
Kepler magnitude is brighter than K p ~ 13.0. The sec- 
ond criterion is adopted since hotter stars are likely to 
have larger V sin I s and so the relative statistic and sys- 
tematic errors in I s become smaller. 

In order to derive spectroscopic parameters for those 
KOI systems, we conducted spectroscopic observations 
with the High Dispersion Spectrograph (HDS) installed 
on the Subaru telescope, located in Hawai'i. The ob- 
servations were performed on June 12, July 23, 24, and 
October 19, 2011 (UT). We employed the I2a observing 
mode, which covers the wavelength regions of 4950 A to 
6150 A (blue CCD) and 6400 A to 7550 A (red CCD). 
The spectral resolutions were set as R ~ 90, 000 on June 
12 and July 24, R ~ 60, 000 on July 23, and R ~ 110, 000 
on October 19, respectively. The resolution for each tar- 
get is summarized in the middle column of Table [T] The 
seeing on each observing nights was typically 0.5 — 0.8 
arcseconds. 

On each of those observing nights, we obtained the ref- 
erence spectrum of the flat lamp transmitted through the 
Iodine cell, using the same setup (spectral resolution). 
These reference spectra were used to reproduce the in- 
strumental profile for each spectral resolution, and then 
play a crucial role in estimating the projected rotational 
velocity V sml s for slowly rotating stars. 

We reduced the raw data of each spectrum using the 
standard IRAF procedure. The resulting spectra have 
the typical signal-to-noise ratio (SNR) of 80 — 120 per 
pixel after extracting the ID spectra. 
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4. ANALYSES AND RESULTS 
4.1. Estimation of Rotational Periods 

In order to determine the periods of photometric vari- 
ations, we used all public data available from the MAST 
archive for the 15 targets. In most cases, quarters 
through 6 were available, corresponding to 500 days of 
observations in total. In this analysis, only the long ca- 
dence observations were used, which proved to be suffi- 
cient to study the flux variability with a timescale of a 
few days. Since the Pre-search Data Conditioning (PDC) 
pipeline is known to remove partially or totally the stel- 
lar flux variability on the timescales we are interested in 
(| Jenkins et al.ll2010bf ). we decided to use the raw flux, 
named SAP_FLUX in the version 2.0 of the FITS files 
delivered by the Kepler team. Note that the newest im- 
provements on th e Kepler pipeline (jStumpe et al.l [20121 : 
iSmith et al.ll2012l ) also seem to allow for a fast identifica- 
tion of the real astrophysical noise, but we simply apply 
the method below in the present paper. 

We used the linear ephemeris of all planet can- 
didates published in the February 2011 data release 
(|Borucki et al.l 120111) to locate all transits and remove 
them from the flux series. Then, using a smoothed ver- 
sion of the flux where each point is the average of the 
previous and next 10 points, we estimate the standard 
deviation of the flux, and we applied 3 sigma clipping 
to the data in order to remove outliers. Raw Kepler 
data suffer from sever al well known instrumental artifacts 
(jJenkins et al.ll2010cD that would not be corrected by the 
3 sigma clipping. Some of them only affect a certain in- 
terval of the observations, for example the changes in 
flux due to temperature drifts after the telescope points 
to the Earth. 

There are also long term trends due to the constant 
movement of the targets on the CCD. In order to remove 
them, we detrend the data using the cotrending basis 
vectors, following the instructions from the Data Release 
12 Notes. These cotrending vectors are constructed to 
capture most of the flux variability caused by instrumen- 
tal artifacts in each CCD. The mean-subtracted, mean- 
normalized flux of each target can be considered as a su- 
perposition of the astrophysical flux of the star plus a lin- 
ear combination of these orthonormal basis of cotrending 
vectors. Assuming that the astrophysical flux is orthog- 
onal to the instrumental noise, we can estimate the co- 
efficients of the linear combination by simply taking the 
vector product of the mean-subtracted, mean-normalized 
raw flux with as many cotrending vectors as needed to 
clean the lightcurve. After this removal, we added back 
the previously subtracted mean, to then normalize by 
the mean each quarter to avoid differences in flux from 
quarter to quarter. 

We applied the above procedure to the 15 targets, and 
studied the final flux lightcurve. Visual inspection of the 
light curves confirmed the existence of flux variability 
due to starspots on 13 out of the 15 cases, with peri- 
odic signals with amplitudes and phases that evolve with 
time. KOI-258 and KOI-302 show a very clear but co- 
herent variability that does not change shape with time. 
Such strictly periodic signals arc unlikely to be caused by 
spots, especially in the case of KOI-258 where the period 
of the signal is equal to the orbital period of the candi- 
date. This strongly suggests that the periodic dimming 



TABLE 1 

Spectral resolutions adopted in spectroscopic 
observations by sljbaru/hds and rotational periods 
estimated by the Kepler photometry - . 



System 


R 


P s (days) 


KOI-4 


~ 110000 


5.65 ± 0.03 


KOI-42 


~ 90000 


20.84 ± 0.37 


KOI-100 


~ 90000 


1.132 ± 0.002 


KOI-257 


~ 90000 


7.846 ± 0.052 


KOI-258 


~ 90000 


too coherent 


KOI-261 


~ 110000 


15.38 ±0.30 


KOI-262 


~ 90000 


8.171 ± 1.218 


KOI-269 


~ 90000 


5.351 ± 0.136 


KOI-279 


~ 90000 


20.92 ± 0.93 


KOI-280 


~ 90000 


15.78 ±2.12 


KOI-302 


~ 60000 


too coherent 


KOI-367 


~ 90000 


27.65 ± 3.56 


KOI-974 


~ 110000 


10.83 ±0.12 


KOI- 1020 


~ 90000 


10.91 ± 1.06 


KOI- 1463 


~ 60000 


6.042 ± 0.042 



of KOI-258 is likely to be caused by a background binary. 

In order to estimate the periods of rotation for the 
rest of the stars, we use a lomb-scarglc algorithm and 
analyze the power spectra of the stellar fluxes. A high 
peak is expected to happen at the period of rotation of 
the star, although a strong peak can also appear at half 
the period if several spots are present. The peak can also 
be rather wide if differential rotation is present, or even 
be composed of several peaks. However, visual inspection 
can also reveal important information about the periodic 
behavior of the flux series, and it helped to identify a few 
cases where indeed the highest peak corresponded to half 
the period of rotation. In these cases, two similar spots, 
with opposite stellar longitudes, induce flux variations 
with twice the frequency. The evolution of the size of 
both spots, which creates asymmetries between the flux 
variations induced by both, allows us to identify them 
unambiguously. 

The final values of the period of rotation and their 
errors are obtained by studying the proximity of the 
strongest peak, and summarized in the right column 
of Table [TJ We adopt the full-width at half maximum 
(FWHM) of the peak as the 1-sigma error, with the cen- 
ter of the interval being the rotation period. Note that if 
a differential rotation produces several peaks with sim- 
ilar amplitudes at similar periods, the resulting merged 
peak in the periodogram is broadened, and the above as- 
signed error becomes large. Thus, the presence of stellar 
differential rotation makes the estimate of the period less 
accurate. 

4.2. Estimation of Spectroscopic Parameters 

Following iTakeda et all (|2002l 120051) . we analyze each 
spectrum and estimate the effective temperature T e g, 
surface gravity log g, microturbulence £, and metallicity 
[Fe/H] by measuring the equivalent widths of Fe I and 
Fe II lines. In order to accurately estimate V sin/ s with 
avoiding any systematic effect, we numerically integrate 
each component on the stellar disk Dopplcr-shifted due 
to stellar rotation and macroturbulcncc. In doing so, we 
also convolve the intrinsic (thermal motion + microtur- 
bulence) profile with the rotational and macroturbulcncc 
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broadening function and the instrumental profile to re- 
produce the observed spectrum. 

For macroturbul ence, we adopt the radial-tangential 
model ; ]( !ravlE(l ()") :. For the given macroturbulence dis- 
persion £rt and the stellar limb-darkening parameter 
(e = 0.6), we determine the best-fit solution of Vsinl s 
for each of the spectra. Since the macroturbulence dis- 
persion £rt is not well understood, especially for hotter 
stars, we try several different values of £rt- For the fi- 
nal ^j^sjah J _jve_ajicj3t the fo llowing empirical expression 
bv lValenti fc Fischer! (|2005l ): 



RT 



3.98 



T cff - 5770 K 
650 K 



km , 



(3) 



and estimate the systematic uncertainty for V sin / s by 
changing £rt by ±15% from Equation §3§ for cool stars 
(Tcff < 6100 K) bas ed on the observed distr ibution of 
Crt (see Figure 3 in IValenti fc Fischer! [20051 ) . For the 
case of hotter stars (T e ff > 6100 K), however, the macro- 
turbulence is not intensely investigated and thus we con- 
servatively estimate the systematic error for l^sin/ s by 
changing £rt by ±25%. The statistical errors in fitting 
each spectrum are generally smaller than the systematic 
errors arising from different values of £rt • 

As for the instrumental profile, we basically adopt 
Gaussian broadening functions whose FWHMs corre- 
spond to each of the spectral resolutions adopted in the 
observations (R ~ 60,000, 90,000, and 110,000). In the 
case of slowly rotating stars (Vsml s < 5 km s _1 ), how- 
ever, we found that the shape of the instrumental pro- 
file, which is slightly different from Gaussian, sometimes 
affects the estimate of V^sin/s. Thus, we convolve the 
actual shapes of the instrumental profile estimated by 
the reference transmission spectrum of the Iodine cell 
fo r those sl ow ro tators (i.e., KOI-261 and KOI-367). 

iTakedal (|1995l ) applied this procedure to the high- 
resolution solar flux spectrum and obtained the solar 



spin velocity of Vsinl s = 2.00 ± 0.34 km s 



Since 



the angular velocity around the solar equator is about 
14° d ay -1 from the observ ations of spatially resolved 
spots (iRuzdiak et al.l I2005T) . we obtain V cq « 2.0 km 



s for the Sun. We also applied this technique to the 
spectrum of HAT-P-11, whose rotational velocity is re- 
ported to be very small (Vsinl s = l.OO+o'gg km s -1 , 
W inn et all l2010bD from the measurement of the RM 
effect, independently of the spectral line analysis. Al- 
though the resolution of HAT-P-ll's spectrum that we 
analyzed is relatively low (R ~ 45000), we obtained 
Vs'ml s = 1.79 ± 0.65 km s _1 from the line analysis 
(including the convolution of the instrumental profile), 
which is consistent with the result estimated by the RM 
measurement. These two tests (for the Sun and HAT-P- 
11) validate the present technique to estimate V sml s for 
slow rotators. Of course, they do not guarantee the ex- 
trapolation of even smaller V sin/ s (< 1.0 km s^ 1 , e.g., 
KOI-261), but for such stars, we can safely rule out a 
large Vs'ml s (> 2 km s" 1 ) with a high confidence level. 

Once T e g, logg, and [Fc/H] arc given, we can estimate 
the stellar age, mass M s , and radius R s for each system. 
We here employ the Yonsei- Yale (Y 2 ) isoc hrone model to 
estimate these parameters (|Yi et al.ll2001f ). The result is 
summarized in Table [2] A visual inspection of the spec- 



FlG. 3. — The slit viewer image for KOI-42. A companion star is 
seen about 2.0" to the north-east of the main star. The companion 
star is on the slit. 




Fig. 4. — The slit viewer image for KOI-279. A companion star 
is seen about 1.5" to the west of the main star. The light from the 
companion is partly on the slit. 



trum indicates that KOI-1020 is a spectroscopic binary, 
which makes it difficult to derive spectroscopic parame- 
ters. In addition, the rotational velocities of KOI-4, 100, 
258, and 366, are so large (Vsm/^ > 30 km s _1 ) that 
their spectra look very flat, which prohibits any reliable 
estimates of the atmospheric parameters (T e s, logg, £, 
and [Fc/H]) with our current SNR's. The best- fit values 
of Vsinls for those systems are estimated as 39.0 km 
s _1 , 32.8 km s _1 , 134.5 km s _1 , and 34.1 km s _1 , for 
KOI-4, 100, 258, and 366, respectively. For such rapid 
rotators, asymmetries in the transit lightcurve may be 
used to deter mine the para meters only if the spin-orbit 
angle is large (Barnes 2009). Such an analysis is beyond 
the scope of this paper, and we do not perform it here. 

Regarding KOI-42 and KOI-279, the estimated values 
of loggr seem too lar ge for stars with T c g of ~ 6500 K 
(e.g., see Figure 2 in iFuhrmannl 1 1 998h . These large val- 
ues of log g may explain the unusually young estimates 
for system ages listed in Tabled This tends to be caused 
by blending of light by a companion star or another back- 
ground source. Thus, we checked the slit viewer images 
for those targets taken simultaneously with the spectra, 
and consequently we found stellar companions for both 
of the systems and those companions seem to be on the 
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TABLE 2 

Spectroscopic Parameters. Starred systems have companion stars, most likely causing contaminations in spectra. 



System T eff (K) 



log 9 



[Fe/H] age (Gyr) M s (M Q ) R s {Rq) V sin h (km s' 1 ) V cq (km s' 1 ) 



KOI-42* 
KOI-257 
KOI-261 
KOI-262 
KOI-269 
KOI-279* 
KOI-280 
KOI-302 
KOI-367 
KOI-974 
KOI- 1463 



6512 ±58 
6218 ±28 
5708 ± 13 
6150 ±53 
6371 ± 50 
6531 ±58 
6047 ± 40 
6616 ± 63 
5667 ± 20 
6385 ± 30 
6578 ± 70 



4.542 ±0.090 
4.286 ±0.055 
4.329 ± 0.030 
3.994 ±0.100 
4.160 ±0.090 
4.425 ±0.100 
4.262 ±0.080 
3.882 ±0.100 
4.279 ± 0.040 
4.058 ±0.050 
3.886 ±0.105 



0.114 ±0.050 
0.154 ±0.031 
0.048 ±0.019 

-0.104 ±0.047 
0.011 ±0.049 
0.313 ±0.063 

-0.258 ±0.034 
0.097 ±0.054 
0.151 ±0.031 

-0.026 ± 0.031 
0.018 ±0.069 



< 0.06 

Z - UU -0.65 
6 52+ - 42 
o Q7+0.60 
-0.52 
1+0.21 
-0.30 



2.21 

< 0.1 
5.16 
1.33 



+0.69 
0.90 
+0.15 
0.13 



6 95+ 77 
"■ 3J _0.79 
9 9 7 +0.12 
-0.14, 
1+0.12 
0.17 



1.48" 1 



1.373 ±0.022 
1 285+ 031 

i - uo '_0.010 
. o 74 +0.091 

1 QO9 + 0.080 

i.ooz_ Q62 
1 450+ ' 022 

1 095+ ' 034 
±.u»o_ Q24 

+0.110 
x - ' lo -0.117 

1 071+ 016 

±,u ' -0.013 
1 4 6 4+0.05l 

-, 79C+O.122 
x - ' zo -0.116 



i 900+11.04') 
i.^ao_ Q Q44 

1 o 47 +0.i05 
1,011 -0 .093 

i - 1DO -0 .046 
1 qro+0.280 
1.3OO_ 25? 

1 616+ ' 227 
±.o±o_ 193 

1 279+ ' 141 
1 ' z,a -0.123 

2.524+^98 
— 0.350 
1 9/1Q+0.055 
l.Z*3_ _p 5g 

1 871+ - 145 

±, °' 1 -0.133 
o 47Q +0.416 
z - 1la -0.359 



13.40 ±0.22 
7.09 ± 0.49 
62+ 109 
10.58 ±0.22 
11.62 ±0.22 
12.14 ±0.27 
3.52 ± 0.50 
16.66 ±0.16 
1.04 ±0.74 
7.13 ±0.49 
9.66 ± 0.36 



3.14 ±0.12 



8.69 



-O.liS 



-0.60 

n on + 0.16 

°- o '- , _0.17 

12.16t 2 ^ 
3.08l°;J? 

,+0.79 
-0.62 



4.12] 

N/A 
2.271 
8.74 4 



7+0.36 
-0.28 
h0.69 
-0.63 

20 76+ 3 ' 49 
zu - '"-3.00 



slit of Subaru/HDS, causing significant contaminations 
in the spectra of KOI-42 and KOI-279 (see Figure |3] 
and 2]) Kepler's photometric aperture is relatively large, 
which means that these companion stars are a big source 
of contamination. This increases the chance that the ro- 
tational periods estimated by spots are biased and might 
reflect the companions' rotational periods. Therefore, we 
simply ignore these two systems in the following analy- 
sis and discuss the stellar inclinations for the rest of the 
systems. 

We also checked the slit viewer images for the other 
targets to see if any contamination sources are located 
around the main objects. Consequently, we found that 
KOI-258 seems to have a companion star, located ~ 1" 
to the east of the main object. This is consistent with 
our expectation from the photometric analysis in Sec- 
tion |4T] In addition, the point spread function (PSF) of 
KOI-1020 looked distorted, which suggests existence of a 
companion or background sources. We could not locate 
any companion (within ~ 1") nor anomalous PSF for the 
other systems. 

We compare our result for the stellar parameters in Ta- 
ble[2]with the public KIC parameters, which are based on 
the photometric analyses. We find that the root-mean- 
square (RMS) differences between our result and KIC 
values for T g and log g are 217 K and 0.231 dex, respec- 
tively. These values are reasonably in good agreement 
with the reported uncertainties for t he KIC parameter s 
(~ 200 K and - 0.4 dex, respectively. iBrown et aL1l2011[ ). 
Moreover, when we remove the possibly contaminated 
systems (KOI-42, 279, and 1463), the RMS differences 
significantly improve and become 101 K and 0.134 dex 
for T e ff and logg, respectively. 

4.3. Evidence of Possible Spin-Orbit Misalignments 

After excluding the systems for which we could not es- 
timate either the rotational period or the stellar radius, 
we are left with KOI-42, 257, 261, 262, 269, 279, 280, 
367, 974, 1463. KOI-42 and KOI-279 were excluded as 
well because of the contaminations by companion stars, 
as mentioned above. In the right column of Table [2 
we also show the rotational velocity at the stellar equa- 
tor V eq based on P s obtained by the Kepler photometry 
and spectroscopically measured R s . In Figure we plot 
the projected rotational velocity V sml s as a function of 
V cq . In this plot, we also show the three lines indicat- 
ing I s = 90°, I s = 45°, and I s = 30°. While KOI-262 
and KOI-280 are consistent with I s = 90°, the other sys- 




FlG. 5. — The estimated V eq and V sin I s . The solid line indicates 
the case that our linc-of-sight is vertical to the stellar spin axis. We 
also plot the two inclined cases in which I B = 45° and 30° . 



terns are inconsistent with I s = 90° within la. In par- 
ticular, KOI-261 and KOI-1463 have significantly small 
stellar inclinations. We note that KOI-262 is a multiple 
transiting system (candidate), which makes the system 
a very important sample to discuss planetary migrations 
(jSanchis-Oieda et al.|[2012]) . 

As for the KOI-1463 system, for which we inferred a 
small sin/,,, the effective temperature of the host star 
(> 6500 K) seems too high to have dark starspots, as 
mentioned in Section [2] This means that the periodic 
flux variation of KOI-1463 may be spurious, implying 
that the flux variations do not originate from the main 
star (with T e s > 6500 K), although it is still possible that 
smaller spots (or spots with a higher brightness) induced 
the flux variation. Indeed the planet-to-star size ratio of 
KOI-1463 is reported to be R p /R s = 0.13655 from the 
Kepler transit lightcurve, and the radius of the transit- 
ing companion becomes 36 .93 jlg'^i^© when we simply 
substitute the stellar radius estimated via spectroscopy. 
An object with such a huge radius corresponds to a very 
late-type star rather than a planet. If this is the case, 
the periodic flux variation may come from the late-type 
companion star, KOI-1463. 01, which may well be active 
enough to have starspots. Although the huge discrep- 
ancy between the public KIC parameter, which reports 
R s = 1.17i?©, and our estimate further supports this 
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scenario, further lightcurve analyses and/or a high res- 
olution imaging are required in order to understand the 
origin of the flux variation of KOI- 1463. 

While I s is just an angle between the stellar spin axis 
and the line-of-sight, a small value of I a (or sin/ s ) for 
transiting systems implies a spin-orbit misalignment, but 
not vice versa. To see this more clearly, let us consider 
the 3D spin-orbit angle -0 using Equation (|2|) again. Re- 
call first that the orbital inclination I Q is expressed as 



cosL, 



l 



e sin w 



1 



(4) 



where 6, a p , e, w are the transit impact parameter, 
planet's semi-major axis, orbital eccentricity, and lon- 
gitude of the pcriastron. Except for extremely eccentric 
planets (e > 0.9), we expect that the right-hand-side of 
Equation Q, apart from the factor b, is of the order 
of R s /a p and b varies from zero to unity by definition. 
Therefore, 



C0Si o < . 



(5) 



Substituting the above relation into Equation ([2]), we 
obtain 



cos 4> = sin I s cos A sin I + cos I s cos I a 

<r ■ t t ^ s t 
5, smi, A cosi«. 



(6) 



Equation (jSJ) gives a lower limit of tp from the observed 
value of I s . For instance, in the case of KOI-261, we 
obtain ip > 61° based on Equation (|6|). 

On the other hand, I s = 90° (i.e., s'ml s = 1) within 
its uncertainty does not necessarily mean a spin-orbit 
alignment. This is analogous to the case of A = 0° for an 
RM measurement, which does not always imply a spin- 
orbit alignment. Statistical treatments are important in 
both cases in order to compare observed distributions 
with planetary migration theories. 



5. DISCUSSION 



5.1. 



Correlation between Stellar Inclinations and Other 
System Parameters 

In the previous section, we have shown that at least one 
system (KOI-261) out of our 15 targets indeed indicates a 
possible spin-orbit misalignment along the line-of-sight. 
KOI-1463 also indicates a small inclination, but given the 
facts that its companion is more likely a star rather than 
a planet and its effective temperature is too high for KOI- 
1463 to have starspots, the inferred rotational velocity 
at the stellar equator V eq is suspicious. Therefore, we do 
not further consider the KOI-1463 system with currently 
available data. Other systems (KOI-257, 269, 974) also 
show possible spin-orbit misalignments with > 2a, but 
those moderate spin-orbit misalignments may be caused 
by stellar differential rotations as we will show later in 
Section I5T31 Before discussing the impact of differential 
rotation, we here discuss the dependences of I s (that tells 
us about spin-orbit relations) on the system parameters. 

Figurc[B]and[7]plot the stellar inclinations sin/ s against 
the stellar effective temperatures T c g and their ages, re- 
spectively. As mentioned in Section [TJ host stars' ef- 
fective temperatures and ages are reported to have sig- 
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Fig. 6. — Correlation between the stellar inclination and stellar 
effective temperature. 



1.2 



0.8 



c 0.6 



0.4 



0.2 



974 




269 



262 



280 




012345678 
age [Gyr] 

Fig. 7. — Correlations between the stellar inclination and stellar 
age. 



nificant correlations with stellar obliquities ([Winn et al.l 
I2010al [TTiaTrcllMTI ). Amongst all, KOI-2 61 does not 
seem to follow the possib le patterns by I Winn et al.1 
(|2010al) and iTriaudl (|201lD . who suggested that spin- 
orbit misalignments are seen around hot, young stars, 
although KOI-261 has only a Neptune-sized planetary 
candidate as shown later in this subsection. 

Next, in order to get a clearer insight into each of the 
systems, we would like to provide a rough estimate for 
planetary parameters by using the public photometric 
data. First, we focus on the semi-major axes of the plan- 
ets. The distance between the planet and its host star 
plays a critical role in discussing planetary migrations, 
tidal interactions, and the resultant stellar obliquities 
with respect to planetary orbits. As usual, the planet 
semi-major axis a p is written as 



a p = {G(M S + M p )} 



1/3 



2/3 



(7) 



where G is the gravitational constant, and P and M p are 
the orbital period and planet mass, respectively. Now 
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TABLE 3 

Correlation between smi" s and Planetary Parameters. 



Planetary Candidate 



sin/. 



KOI-257.01 
KOI-261.01 
KOI-262.01 
KOI-262.02 
KOI-269.01 
KOI-280.01 
KOI-367.01 
KOI-974.01 
KOI-1463.01 



U.814_ Q7 g 

q 162 +0.283 



P (days) (adopted) R p /R s (adopted) 



a P (AU) 



0.869 
0.869 
0.760 
0.846 
0.448 
0.813 
0.465 



0.162 
+0.192 
0.166 
+0.192 
0.166 
+0.105 
0.097 
+0.201 
0.175 
+0.336 
0.321 
+0.085 
0.079 
+0.080 
0.069 



6.883403 ± 0.000012 
16.238480 ± 0.000019 
7.8125124 ± 0.000052 
9.376137 ± 0.000056 

18.01134 ± 0.00022 
11.872914 ± 0.000023 
31.578680 ± 0.000018 

53.50607 ± 0.00061 
N/A 



0.02052 ±0.00015 
0.02431 ± 0.00033 
0.01074 ± 0.00015 
0.01362 ± 0.00030 
0.01074 ± 0.00019 
0.01972 ± 0.00073 

0.0420 ± 0.0038 
0.01353 ± 0.00014 
0.13655 



0769 +() ' ,,UUV 
u - u ' Da -0.0004 

0.1278 ± 0.0004 

o.ossetnSi? 

0.0967±°:gS|S 



0.1497"* 



7 +0.0029 
'-0.0022 
1050+ - 0010 

u.iuou_ 0008 

o 9nno+ - 0010 
u.zuuu_ 0008 

SI 55+0 0037 
U.OlOO_ 0034 

N/A 



°- ul +0.21 

-0.33 

'+0.30 

2 91 - 42 
z ' M± +0.38 



2.29] 



1.89] 



-0.27 
'+0.23 
o ^c — 0.32 
A - '°+0.28 
c cq-0.59 
°' DM +0.57 
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e 0.6 




0.01 0.1 

semi-major axis [AU] 

FlG. 8. — Correlation between the stellar inclination and semi- 
major axes of planets. 
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Fig. 9. — Correlation between the stellar inclination and plane- 
tary radii. 



that we have a better estimate for the stellar mass M s 
from spectroscopy, we substitute M s in Tabled] into the 
above equation. Adopting the public orbital period for 
each of the KOI planetary candidates from the MAST 
archive (as of March 2012), and neglecting the term of 
M p /M s (< 0.001), we obtain the semi-major axis a p 



shown in Table [3] We plot the stellar inclination sin I s for 
each planetary candidate as a function of the semi-major 
axis a p in Figure [8] 

Finally, we deal with the planet radius. From the tran- 
sit lightcurve, it is generally possible to estimate the size 
of the transiting planet relative to the stellar radius, but 
since the systems on which we are focusing are only can- 
didates, we need to keep in mind that lightcurve may 
be blended with the flux from a binary companion or 
other background sources. In particular, because of the 
relatively large aperture size of the Kepler photometry, 
some of the systems may well be blended. To quantify 
the rejection level of blending scenarios, we estimated 
the upper limits of the secondary peaks in the spectral 
lines. Assuming the same line profile, we tried to fit a 
secondary peak for each of the observed spectra with a 
differing flux ratio and Dopplcr shift. As a result, we typ- 
ically obtained rejection levels (ler upper limit) of 0.01 
to 0.11 for the flux ratio between the two peaks. The 
rejection levels are worse for rapidly rotating stars (i.e., 
KOI-262, 269, 974, and 1463) because of the shallow line 
profiles. Note that our spectra typically have an S/N of 
~ 100, so that the determinations of the upper limit are 
dominated by the photon noise. Further observations are 
required in order to rule out the binary scenario. Here, 
we simply adopt the planet-to-star size ratios R p /R s re- 
ported by the Kepler team, and infer the planet radii 
using the spectroscopically estimated stellar radii. 

The result is also shown in Table [3] and its correlation 
with sin I s is plotted in Figure |H1 Most of the candi- 
dates seem to be super-Earths or Neptune-sized plan- 
ets. As for the two Neptune-sized planets orbiting rel- 
atively cool stars, the orbit of KOI-261.01 seems to be 
inclined with respect to the spin axis of the host star, as 
in the case of the super-Neptune HAT-P-llb, for which 
a significant spin-orbit misalignment aroun d a cool star 
was found through the RM measurements (|Winn et al.l 
I2010bt IHirano et all [2CJlTah and a precise photometric 
analysis of the effect of spots in the transit light curves 
(jSanch is-Oicd a fc Winnll201lD . As we have noted in Sec- 
tion [TJ giant planets (including Neptunc-sized ones) dis- 
covered inside of the snow line (~ a few AU) should have 
experienced planetary migrations. The possible spin- 
orbit misalignments in the KOI-261 system suggests that 
in some of the systems with close-in Neptune-sized plan- 
ets, planet-planet scatterings or other dynamical pro- 
cesses such as the Kozai cycles may have played impor- 
tant roles during their formations and evolutions. 



Measurements of Stellar Inclinations for KOI Systems 



9 



TABLE 4 

Rotational periods estimated by the empirical relation. 



System 


P 3 (days) 


P 3 , model (days) 


Flux Variability (%) 


KOI-257 


7.846 ± 0.052 


6 56 +tJ - 5y 

D - oo -0 .79 


0.1511 


KOI-261 


15.38 ±0.30 


25.98±°'^ 


0.3357 


KOI-262 


8.171 ± 1.218 


6.5ll?;|° 


0.0520 


KOI-269 


5.351 ±0.136 


5.061°;^ 


0.0213 


KOI-280 


15.78 ±2.12 


20.33l24 7 


0.0318 


KOI-367 


27.65 ±3.56 


25.8OI208 


0.1161 


KOI-974 


10.83 ±0.12 


4 26+ ' 56 


0.0355 


KOI- 1463 


6.042 ± 0.042 


1 64+ 1 ' 36 


0.0378 



1.002 



1.000 



£ 0.998 - 



0.996 




5.2. Comparison with Empirical Estimates for P s 

We have estimated the stellar rotational period P s 
based on the periodic analysis of the Kepler photome- 
try While the current method is valid as long as the 
periodic flux variation comes from the spot on the stellar 
surface, the variation may be due to some other sources 
such as flux variations of background stars. In order to 
check the reliability of P s that we estimated, we apply the 
Kyroch ronological method as employed by iSchlaufmanl 
(|2010t) . We simply adopt the same empirical relation for 
the rotational period as a function of the stellar mass M s 
and age. Substituting M s and ages lis t ed in Table [2] into 
the empirical formula by ISchlaufmanl (|2010T l , we obtain 
the modeled rotational periods for the KOI systems in 
Table [31 for which secure spectroscopic parameters are 
obtained. The derived periods -P s ,modci are listed in Ta- 
ble H] along with the observed rotational periods P s . The 
uncertainty in P s . mo dci comes solely from the uncertain- 
ties in M s and age inferred from spectroscopy, and does 
not include any systematics in the empirical relation. 

The comparison between P s and P s , model indicates that 
for most of the systems the rotational periods derived 
from stellar spots are in reasonably good agreement with 
the empirical ones, while KOI-261, 974, 1463 show some 
discrepancies. This suggests that the estimations of ci- 
ther of P s or spectroscopic parameters may be wrong for 
those systems. It should al so be emphasized, however, 
that the empirical model by Schlaufman (2010) is rather 
simplified, and it is difficult to evaluate the systematic 
errors caused by adopting that empirical relation. For in- 
stance, we note that KOI-1463's host star is so massive 
that the empirical model may have significant systemat- 
ics because of t he lack of the samp led stars in that region 
(see Figure 2 in ISchlaufmanl l20T0l) . 

There is a non-zero probability that the flux variabil- 
ity is caused by a background star, a probability that 
is higher, the lower the observed variability is. To in- 
vestigate this, we estimate the flux variability by taking 
the detrended mean normalized flux, and eliminating the 
5% highest and the 5% lowest values. This step helps to 
remove the effect of outliers or artifacts that might still 
remain after the detrending process. The flux variability 
is then defined as the range of values of the flux. In Table 
IU we show the values of the variability for each system. 
As one can see, KOI-974 and KOI- 1463 show a very low 
level of variability, with a higher false positive proba- 
bility, whereas KOI-261 is the most active star, which 
shows that the period of rotation has been calculated 
more robustly. Note that assuming that the flux vari- 



420 440 460 480 500 

BJD-2454900 [days] 

Fig. 10. — The black dots represent the corrected flux series of 
KOI-261 on a segment of 90 days. In this segment, the star has a 
peak-to-pcak variability of up to 0.6%. The blue and red arrows 
point to the flux minima generated by two different active regions. 
The active region represented by blue arrows seems to disappear 
after three rotation periods, whereas the one represented by red 
arrows seems to reach its maximum size in the middle of the ob- 
servations. In both cases, the flux minima recur with a periodicity 
of 15 days, confirming the value obtained with the Lomb-Scarglc 
pcriodogram. 

ation in KOI- 1463 comes from the low-mass companion 
star and that the flux ratio between them is 50 to 100, a 
0.0368% variability in the total flux corresponds to 1.89% 
to 3.78% variability in the companion's flux, which seems 
to be reasonably caused by starspots on an active late- 
type star. 

Since KOI-261 is apparently the most important sys- 
tem in our sample that most likely to have a spin-orbit 
misalignment, we pay a special attention to this system. 
If we adopt the modeled rotational period for KOI-261 
as shown in Table 0] we obtain V^ q = 2.27 ±0.10 km 
s _1 , which agrees with its projected rotational velocity 
of V sin/ s = 0.621q km s" 1 within 2<r. However, as we 
have shown in Figure [2] the pcriodogram clearly shows 
the single strongest peak around 15.5 days. The absence 
of a peak around 31 days implies that the rotational pe- 
riod is securely derived and does not reflect half the pe- 
riod. We also show in Figure [TU] that there exist two 
active regions on the stellar surface manifested in KOI- 
261 's lightcurve. The reason for the disagreement be- 
tween the observed and modeled rotational periods for 
KOI-261 is unknown, but it is general believed that such 
active stars as KOI-261 rotate faster. Future confirma- 
tion and characterizations of this system is particularly 
intriguing. 

5.3. Impact of Differential Rotation 

So far, we have discussed the stellar inclinations / s 
assuming no differential rotation of the planet hosting 
star. In reality, however, stars may have differential ro- 
tations, which are supposed to be related to the origin 
of starspots, and the rotational velocity at the equator 
may be different from that estimated at the location of 
the star spot. Here, we discuss the impact of stellar dif- 
ferential rot ations on our estimate of I s . 

Following iReiners fc Schmi tt (2003), we model the an- 
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Fig. 11. — The simulated result for the line profile fitting in 
the presence of differential rotation. Wc plot the best-fit values 
for (V sin/ s ) spcc obtained by fitting the line profiles of differential 
rotation with the rotational kernel for a "rigid rotation" , for various 
input stellar inclinations (sin/ s )i n . In this plot, we employ V eq = 
10 km s — 1 and a = 0.2, and the black line represents the expected 
V sin / s for rigid rotation. The best-fit results for (Vsin/ S ) sp ec are 
deviated from the cases for rigid rotation by ~ 10%. 

gular velocity f2 as a function of the latitude I on the 
stellar surface as 

ft(l) = fi eq (l -a sin 2 /), (8) 

where £l eq is the angular velocity of the star at the equa- 
tor. The degree of differential rotation a is about 0.2 
for the case of the Sun. In the presence of a differential 
rotation, it is more explicit to rewrite Equation ([I} as 

sin7 3 = (V** 1 -)*- (9 ) 
i? s £l oq (l — a sin I) 

Since we do not have information on I where the stellar 
spot that we observed is located, the uncertainty for I 
results in a systematic error for sin I s . A fortunate situa- 
tion is, however, the spectroscopically measured V ainl s 
in the numerator of Equation ^ is likely estimated to be 
smaller in the presence of differential rotations (a > 0). 
This is because the solar-type differe ntial rotation leads 
to a "sharper" spectral line profile ()Reiners fc Schmittl 
I2003T ) , which results in a smaller value for V sin I s when 
fitted by the rotational broadening kernel for "rigid ro- 
tation" . Since the denominator of Equation © also be- 
comes smaller in the presence of differential rotation, the 
impact of differential rotation tends to be more or less 
compensated. 

In order to quantitatively evaluate Equation ([§]), we 
perform the following simple numerical simulation. First, 
we generate the mock line profile by convolving a single 
Gaussian function with a rotational broadening kernel 
(including macroturbulence) for a differentially rotating 
star. The input parameters to create one line profile 
are the rotational velocity at the equator V eq = R s £l cq , 
differential rotation parameter a, and input stellar incli- 
nation (sin/s)^. Then, we fit the resultant line profile 
with a convolution function between a Gaussian and a ro- 

8 The other parameters such as the original Gaussian width and 
macroturbulence dispersion are fixed in the simulation. 



Fig. 12. — Same as Figure [5] after corrected for the impact of 
differential rotations. 

tational broadening kernel for a rigid body, and estimate 
the best-fit value for (V sml s ) spec . 

As a result, in the case of a = 0.2, the best-fit values 
for (V sin I s ) spe c are smaller than the input values of V cq ■ 
(sin/ s )i n by ~ 10%, for various cases of (sin/ s ) in (see 
Figure II 1 [) . This result implies that changing sin I from 
to 1 in Equation ((9]) leads to the systematic errors in 
the output sin/ s of ~ ±10%. We check it for various 
values of the input parameters V^q, a, and (sin/ s )i n and 
find that the relative systematic error in estimating sin I s 
from Equation ^ is approximately 0.5a. 

In the case of the Sun, however, the latitudes at which 
we observe the spots are confined to the relatively narrow 
bands of 5° < |Z| < 40°. The active latitude is known 
to move toward the solar equator with a cycle of ~ 11 
years, making the w ell-known "butterfly diagram" (e.g., 
iRuzdjak et aLl[2005T l. 

In order to correct for the impact of differential rota- 
tion, we her e adopt the following e mpirical model for a 
reported bv lCollier Cameron] (|2007f ): 

aQ cq = 0.053(T c ff/5130) 8 - 6 . (10) 

This expression, which is based on observations of 
Doppler imaging or Fourier analyses of the rota- 
tional broadening kernel, claims that differential rota- 
tions are stronger for hotter stars such as o ur targets 
([Collier Cameronll2007t fPunstone et al.ll2008f ). We per- 
form the same simulation using the mock line profile in 
the presence of differential rotation for each of the seven 
systems, and then correct for the impact of differential 
rotation. We evaluate the systematic errors in V eq as- 
suming that the spots are located at |Z| = 20° ± 20° , 
imitating the case of the Sun. 

Figure [T2l plots thus corrected V cq - ysin/ s diagram. 
In contrast to Figure El KOI-257, 269, 367, and 974 are 
consistent with I s = 90° with ~ la, while KOI-261 still 
seems to have a spin-orbit misalignment. This implies 
that differential rotations can partially or totally explain 
the moderate and apparent spin-orbit misalignments for 
KOI-257, 269, 367, and 974, although we do not conclude 
so in the present paper because of the rather rough as- 
sumptions wc adopted for the correction of differential 
rotation. We note that if this is the case, it is equally in- 
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teresting in the sense that the present method can reveal 
the signature of differential rotations of distant transiting 
planetary systems. 

5.4. Comparison with the RM Measurement for 
Kepler-8 

Finally, we briefly discuss the application of the present 
technique to the confirmed Kepler systems. Specifically, 
systems where the RM effect has been measured are par- 
ticularly interesting since we can directly compare be- 
tween the sky-projected and line-of-sight spin-orbit an- 
gles. Among the confirmed Kepler systems with the RM 
measurements, Kcpler-8 shows a periodic flux variation, 
and we obtain the rotational period of P s — 7.5 ± 0.3 
days. This period, along with the stellar radius of 
R s = 1.486^;^i? (jJenkins et all l2010a| ). results in 
V cq = 10.0±0.6 km s -1 , which is in good agreement with 
the proj ected rotational veloc ity Fsin/ S = 10.5 ± 0.7 
km s" 1 (jJenkins et all l2010af l. Although I s ps 90° is 
only a requirement for a spin-orbit alignm ent, this result 
does no t support the RM measurement bv I Jenkins et afl 
(I2010afl (A = -26.4° ± 10.1°). See also lAlbrechti (120121) . 
who suspected the uncertainty of A for Kcpler-8 to be 
underestimated. 

6. SUMMARY AND FUTURE PROSPECTS 

We have investigated the stellar inclinations for KOI 
planetary system candidates on the basis of the detailed 
analyses of photometric variation due to stellar spots and 
spectroscopic measurements of V sin I s . We have found 
that at least one system, KOI-261, exhibits a strong sig- 
nature of a possible spin-orbit misalignment along the 
line-of-sight. The planetary candidate KOI-261. 01 is a 
Neptune-sized one (R p ps 3i?©), with a moderate orbital 
distance (a p ps 0.13 AU). If this system is confirmed and 
eventually turns out to be misaligned, KOI-261. 01 will 
become the smallest planet ever reported to have a spin- 
orbit misalignment, which makes it an important sample 
to test and discuss planetary migrations in an extended 
parameter space. KOI-1463 also shows a small stellar 
inclination, though the size of the transiting companion 
may correspond to that of an M star and the periodic sig- 
nal may come from the companion. The results for the 
other systems are ambiguous but they may be interpreted 
as either mildly misaligned or differentially rotating. 

One of the next tasks is to increase the number of 
samples and further discuss the correlations between I s 
and other system parameters, for each of the Earth-sized 
planet population and giant planet population. This will 
require to choose fainter stars, unless new updates of the 
catalog bring many candidates around variable hot stars. 
Even with the currently available dataset, it is relatively 
easy to plan a more efficient observing run, giving higher 
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priority to stars with a higher variability and ruling out 
those ones with no variability or with signals that are too 
coherent to be caused by spots. 

It is also important to refine the measurements of I s . 
We are potentially able to do so if we combine the spec- 
troscopic measurement with photometric data. For in- 
stance, the ratio of the planet semi-major axis to stellar 
radius is available via transit photometry, but we did not 
use it in estimating the stellar radius R s . This is simply 
because a/R s generally depends on the orbital eccentric- 
ity e, and e cannot be determined by the transit alone 
unless the secondary eclipse is seen in the lightcurve. If 
radial velocities for the systems presented in this paper 
are followed-up by some future observations and their or- 
bital eccentricities are measured, we are able to constrain 
R s more precisely and the estimates for sin/, would be 
significantly improved. 

Finally, we stress that all the systems that we ana- 
lyzed are still candidates for having planetary compan- 
ions, and therefore we must be careful in discussing the 
evolution history of planetary systems with the present 
technique. However, we can safely say that possible mis- 
aligned systems suggested by our analysis would become 
very interesting targets for future RM measurements. 
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